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Abstract 

In any geometrically nonlinear quadratic Cosserat-micropolar extended continuum model for¬ 
mulated in the deformation gradient field F := Vy : fl —)■ GL“''(n) and the microrotation field 
i? : 0 —>■ SO(n), the shear-stretch energy is necessarily of the form 

:= /r ||sym(i?^n - 1)||^ -h /Tc ||skew(i?^F - 1)||^ , 

where /r > 0 is the Lame shear modulus and /ic > 0 is the Cosserat couple modulus. In the 
present contribution, we work towards explicit characterizations of the set of optimal Cosserat 
microrotations argmin^g as a function of F G CL'''(n) and weights /x > 0 

and fic ^ 0. For n > 2, we prove a parameter reduction lemma which reduces the optimality 
problem to two limit cases: = (Ijl) and = (IjO)- In contrast to Crioli’s 

theorem, we derive non-classical minimizers for the parameter range /i > /Xc > 0 in dimension 
n = 2. Currently, optimality results for n > 3 are out of reach for us, but we contribute explicit 
representations for n = 2 which we name rpolar^^^(F) G SO(2) and which arise for n = 3 by 
fixing the rotation axis a priori. Further, we compute the associated reduced energy levels 
and study the non-classical optimal Cosserat rotations rpolar^^^(Fy) for simple planar shear. 
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1 Introduction 


In 1940 Guiseppe Grioli proved a remarkable optimality result [B] for the polar factor in dimension 
n = 3. To state his result, we denote by Rp(-F) € SO(n) the unique orthogonal factor of F G 
GL'*"(n) in the right polar decomposition F = Rp(F) U{F) and by U{F) = Rp(F)^F = y/F'^F G 
PSym(n) the symmetric positive definite Biot stretch tensor. In Grioli proved the special case 
n = 3 of the following theorem: 

Theorem 1.1 (Grioli’s theorem [B1IT2E])- Letn > 2 and ||X||^ := tr the Frobenius norm. 

Then for any F G GL“''(n), it holds 

argmin lll^ = {Rp(F)}, and thus min lll^ = ||[/— 1 ||^ . (1-1) 

flGSO(n) i?GSO(n)" 

The optimality of the polar factor Rp(F) for n = 3 generalizes to any dimension n > 2, see, 
e-g-, [l2], and it is this more general theorem to which we shall refer as Grioli’s theorem in this 
present work. A modern exposition of the original contribution of Grioli has been recently made 
available in m- 

In contrast to Grioli’s theorem, it has been noted in [15] that the polar factor Rp(F) is not 
necessarily optimal for a more general formulation of Grioli’s theorem with weights. Hence, our 
main objective is to make progress on 

Problem 1.2 (Weighted optimality). Let n > 2. Compute the set of optimal rotations 

argmin F) := argmin |/r ||sym(ii^F — 1 ) ||^ + p,c ||skew(i?^F — (1.2) 

flG SO(ra) ’ flG SO(n) ^ ^ 

for given F G GL“''(n) and weights ^ > 0,/ic > 0 such that n ^ /ic- Here, sym(A) := ^(X + X^) 
and skew(X) := ^(X — X^). 


For p, = we recover Grioli’s theorem. This case is well understood. However, no systematic 
analysis of the proposed weighted optimality problem seems to exist in the literature. 


Note that Problem |1.2| is of independent interest in the mechanics of micropolar media and Gosserat 
theory in particular. The weighted strain energy contribution arises in any geometrically 

nonlinear Gosserat-micropolar model proposed in, e.g., m, cf. in m [ini usi nn nn 113 and for 
any 6-parameter geometrically nonlinear Gosserat shell model proposed in [IlllIllMlIlSj- In the 
beforementioned material models determines the shear-stretch energy contribution. To 

obtain a full Gosserat continuum model, needs to be augmented by a curvature energy 

term [^ and a volumetric energy term, see, e.g., [T^or [T5] . However, Problem 1.2 reappears as a 
limit case for vanishing characteristic length Lc = 0 j In this scenario, the curvature term vanishes 


identically and the set of global minimizers for (1.2) admits an interpretation as optimal Gosserat 
rotations for given deformation gradient F G GL~'’(n) and material parameters p, and pc- Grioli’s 
theorem implies that for p = pc the optimal Gosserat rotation is uniquely given by Rp(F), i.e., 
it is the orthogonal part of the right polar decomposition F = Rp { 7 . The corresponding reduced 
energy level is the Biot-energy, see [T5] . 

The suitable choice of the so-called Gosserat couple modulus Pc > 0 for specific materials or 
boundary value problems is an interesting open question |14] . In particular, a strictly positive 
choice /ic > 0 is debatable [T3] and a better understanding of the limit case pc = 0 is hence of 


interest. This further motivates to study the weighted formulation stated in Problem 1.2 


We want to stress that although the term W^^^^ subject to minimization in (1.2) is quadratic in 
the nonsymmetric microstrain tensor t/ — 1 = F^F— 1 (see, e.g., 0), the associated minimization 
problem with respect to R is nonlinear due to the multiplicative coupling of R and F and the 
geometry of SO(n). 


^For this interpretation to work, we have to assume that the volume term decouples from the microrotation R, 


(det[r/] - 1)^ -t ( 


V det[[/] 


- 1 


which seems natural, since it is trivially satisfied in the linear 


e.g., := f 

Gosserat models [20 I14II19| . The linear elastic energy density is of the form A) := ||sym(Vu — A)\\^ + 

fic ||skew(Vu — A)\\'^ + ^ tr [(Vu — A)]‘^. Here, for all admissible values of /r > 0 and fic > 0, the unique minimizer 
is given by A = skew(Vw). 
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Remark 1.3 (Existence of global minimizers). The energy is a polynomial in the 

matrix entries, hence G C°°(SO(n), M). Further, since the Lie group SO(n) is compact and 

dSO(n) = 0, the global extrema are attained at interior points. 


The previous remark hints at a possible solution strategy for Problem |1.2[ Suppose that we 
succeed to compute all the critical points i?crit € SO(n) of Then, if possible, a 

direct comparison of the associated critical energy levels ; F) might allow us to identify 

the energy-minimizing branches for given parameters E,/i and p,c- Note that any minimal branch 
coincides with the reduced Cosserat shear-stretch energy 


IT-" : GL+(n) 




W: 


red 

/1,/ic 


(F):= min 

flGSO(n) 


(1.3) 


At present a solution for the three-dimensional problem (let alone the n-dimensional problem) 
seems out of reach for us. Therefore, we restrict our attention to the planar case where we can 
base our computations on the standard parametrisation 


R : [-TT, tt] ^ SO(2) C R(a) 


f COS a — sin a 
\^sina cos a 


(1.4) 


by a rotation anglej^ 

Regarding the material parameters, we prove that, for any dimension n > 2, it is sufficient to restrict 
our attention to two parameter pairs: (/r,/ic) = (1,1)) the classical case, and {fJ., Hc) = (1,0), the 
non-classical case. We shall see that, somewhat surprisingly, the solutions for arbitrary fi > 0 and 
> 0 can be recovered from these two limiting cases. A large part of this paper is dedicated to 
the discussion of the non-classical choice of material parameters in the planar case. 

Problem 1.4 (The planar minimization problem). Let F G GL’'’(2), fx > 0 and fXc^ Q- The task 
is to compute the set of optimal microrotation angles 


argmin 

Qt G [ —TTjTt] 


(v^ skew). 


f COSO 

l^sin a 


— sina\^ 

(Fn 

cos a J 

[F21 


Fi2'\ 

F 22 J 



(1.5) 


It turns out that there are at most two optimal planar rotations in what we will discern as the 
non-classical parameter range fx > p,c > 0. Both of them coincide with the polar factor Rp(E) 
in the compressive regime of A G GL'''(2), but deviate elsewhere; see Section]^ We denote the 
explicit formulae for the optimal Cosserat rotations (i.e., solutions to ( |i^ ) by rpolar^_^^(F) and 
the respective rotation angles (i.e., solutions to (1.5)) by (possibly multi-valued). The 

computation of global minimizers in dependence of F is not completely obvious even for the reduced 
planar case. We hope that the discovered mechanisms will help to understand the cases n > 3 
eventually. 


This paper is now structured as follows: we first prove a dimension-independent parameter trans¬ 
formation lemma for fx and fXc in Sectionj^which allows us to focus on a classical {fx, /ic) = (1,1) and 
a non-classical limit case {p,,fXc) = (1,0)- In Section]^ we compute the optimal planar Cosserat 
rotations and associated formally reduced energies, first for the classical and non-classical limit 
cases for p and /ic, then for all admissible values of the material parameters. Finally, the optimal 
rotations for the non-classical limit case (/i,p.c) = (1,0) are specialized to the particular case of 
planar simple shear in Section A short appendix provides some elementary but useful matrix 
identities for n = 2. 


2 Parameter reduction for /x and fic 


The Cosserat shear-stretch energy (1.21 is parametrized by two material parameters: ^ > 0 and 
fie ^ 0- It can be shown that fi must coincide with the classical Lame shear modulus of linear 


^Since the boundary of SO(n) is empty the rotation ilcrit is a critical point at F G GL'*"(n) if and only if 
^ ~ every smooth curve of rotations R{t) : (—e,e) —>■ SO{n) passing through /?{0) = 

-^crit • 

^Note that tt and — tt are mapped to the same rotation. In this text, we implicitly choose tt over — tt for the 
rotation angle whenever uniqueness is an issue. 
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elasticity. The interpretation of the Cosserat couple modulus /Tc is, however, less clear. It seems 
that all measurements that can be found in the literature lead to inconsistencies for /ic > 0, see, 
6-g-, |f4j . and so a better understanding of these material parameters and their interaction is of 
interest. In this section, we contribute some helpful representations of the Cosserat shear-stretch 
energy in terms of powers of tr [ [/ ] = tr . 

Let us introduce the following equivalence relation for continuous functions /, g : X —> M dehned 
on a compact set X\ 

f 9 argmin/(a;) = argmin 5 (a;) . (2.1) 

X G X X £ X 

An important example is given by / A/ -I- c for A > 0 and c both independent of a;j^ In the 
present work, we consider minimization w.r.t. X = SO(n) and we shall simply write f g instead 
of / ~SO(n) 9 - 

In this section, we show that it is sufficient to restrict our attention to two representative pairs of 
parameters, the classical limit case {g, Hc) = (1,1) and the non-classical limit case (/i,/ic) = (1,0). 
The solutions for arbitrary admissible fi and fXc can then be recovered from these two limiting 
cases by a suitable transformation, see Section [3)4] for details. To this end, we first introduce the 
following 

Definition 2.1 (Parameter rescaling). Let fj, > Hc > 0. We define the singular radius by 

Pa a '■= —— > 0 , and further define Xa a ■= = —-— , ( 2 . 2 ) 

P-Pc ’ Pi,o P-Pc 

as the induced scaling parameter. Note that pi q = 2 and Ai^o = 1- Further, we define the 

parameter rescaling given by 

■■= A-^A = eGL+(n). (2.3) 

For p > 0 and pc = 0, we obtain = F, i.e., the rescaling is only effective for pc > 0. We can 
now state the key result of this section which is independent of the dimension. 

Lemma 2.2 (Parameter reduction). Let n>2 and let F G GL'^(n), then 

Pc>p>0 W^,^^(i?;F) ~ Wi^i{R-F), and 

P>Pc>0 W^,^^(i?;F) ~ Wi,o(i?;F^,^J . 


We have to defer the proof for a bit, since we need some more preparations. 


In particular, the foregoing Lemma [2^ implies that we can focus on the classical and non-classical 
limit cases. Once these are solved, the solutions for general values of p and Pc can be recovered. 
Note that for the non-classical limit case {p, pc) = (1,0), the Cosserat shear-stretch energy defined 
in (1.21 which is subject to minimization, takes the explicit form 


Wi^g: SO(n) X GL'''(n) 


p + 
^0 : 


Wi,o(i?;F) = ||sym(i?^F-l) 


(2.5) 


Remark 2.3 (Immediate generalizations). The parameter reduction in Lemma 2.2 and the 


optimality results for n = 2 in Section both extend to bijective parameter transformations 
T : GL'''(n) -G GL“''(n). For example, t{F) := Gof (F)^, Cof (A) := det[A] arises in the 
cofactor shear energy 


W;°l{R-F) : SO(n) x GL+(n) ^ K+ , 

:= P II sym(Gof (F^F) - + p, || skew(Gof (F^F) - l)f 

= p II sym(F^T(F) - 1)||2 + p, || skew(F^T(F) - l)f . 


Our results extend to this case in the natural way, i.e., the optimal rotations are (rpolar^ or)(F). 


In what follows, we take a detailed look at the parameters p and p^ and the role they play in the 
weighted minimization in Problem |1.2[ 

"^Note that / and A/ + c share the same critical point structure, whereas / ~ g does not imply this. 
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2.1 Reduction of the classical parameter range: > 0 

For the classical parameter range /Xc > /i > 0, we essentially rediscover Grioli’s theorem (see The¬ 
orem o on the optimality of the polar factor Rp(F"), since the minimization problem reduces to 
the limit case (/x,/ic) = (1, !)■ 


Proof of Lemma 2.2 (first part). For /Xc > /x > 0, we have /x^ 
lower bound 


/X > 0 which gives us the following 


:= /X II sym(i?'^F) - l|p /x^ || skew(R^F - l)|p 
= fj.\\R^F - l|p -I- {pc - p) II skew(i?^F - l)|p 
PpWR^F-^W^ = pWip{R-F). 

Since Wp_p^(Rp(F") \F) = ^fFip(Rp(F) ;F), it follows that ^ IFi^i for the entire classical 

parameter range. ■ 


In passing, we have proved the following immediate generalization to Grioli’s theorem. 

Corollary 2.4. Let Pc> p> 0 and F G GL^(n), then 
argmin = argmin {p || sym(i?'^F’) - l|p -|- II skew(R^F - 1)|P} = {Rp(F’)} , 

iXGSO(n) HeSO(n) 

i.e., the polar faetor Rp(F’) is the unique global minimizer. I 


2.2 Reduction of the non-classical parameter range: fj, > Hc > 0 


A sparking idea which enters the proof of the second part of Lemma 2.2 is due to M. Hofmann- 
Kliemt (then at TU Darmstadt [7]) who contributed to the study of the influence of the parameters 
p and Pc by spotting the applicability of the following elementary identity [8]: 

Lemma 2.5 (Expanding the square). Let R G SO(n) and F G GL^(n), then the following identity 
holds: 


tr 


{R'^F - Pp.p.l) = tr [(R^F)2] - 2pp,^^tr [R^F] + p^,p^tr [1]. 


Proof. 


tr 


{R^F - J = tr [{R^F - {R^F - 

= tr [{R^Ff - 2 R'^F + l]. 


( 2 . 6 ) 


(2.7) 


The claim follows by linearity of the trace operator. 


This leads now to a reduction of the minimization of the energy to Wi^o for the non-classical 

parameter range. The main ingredient is the rescaling of the parameter space GL~'’(n) in Definition 

o 


Proof of Lemma 2.2 (second part). We proceed by successive term expansion, gathering the con- 

( 2 . 8 ) 


tributions which are constant with respect to R at each step. To this end, we split 

II sym{R^F - l)|p + || skew(i?^F - l)|p 


= : I =: II 

and simplify the summands I and II separately. For the first term, we get 

I = M II sym {R^F - 1)||2 = | (||Ef + 2||lf + tr [{R^Ff] - 4tr [R^F]) . 


(2.9) 


Similarly, for the second term 

II = /Xc||skew(i?^F)f = ^ {r'^F-F^R, R^F-F^R) = ^ (||Ff - tr [(R^F)^]) (2.10) 
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is obtained. Summation of I and II while shifting all terms constant in R to the right yields 


; F) = I + II = tr [(i?^F)2] - 2^ tr [R^F] + |lFf + ||lf. (2.11) 

We shall collect all terms which are constant with respect to i? in a sequence of suitable constants, 
starting with := || J'lP +^||l|p. This yields the expression 

= ^i^tr[(i?^F)2] -2/itr[i?^F] +c^l(F). (2.12) 


Introducing the singular radius 
follows 


from Definition 2.1 


we can write the preceding equation as 




^ tr [(i?^F)2] - 2p tr [R^F] + ^(F) , 


(2.13) 


which inspires us to define a rescaled energy 


W^.^,(F;F): SO(u) x GL+(n) ^ K+, W^,^,(F;F) := 




W^.^„(F;F) . 


(2.14) 


We now expand using (2.13) to get 


W^,^^(F;F) = ^ W^.^,(F;F) = ^ f ^tr [(F^F)^] - 2p tr [F^F] +cW^(F)) 

= tr [(F^F)2] - 2p^,^^ tr [F^F] + c^^l^iF) , (2.15) 

with c|t^|i„(F) := c^^)i^(F), and observe that W^^^^ and the rescaled energy share the 

same local and global extrema in SO(n). This gives us W^_^^ ~ which also holds for the 

particular choice of parameters p = 1 and Hc = 0, i.e., Wi^ ^ Wi,o. For the latter specific choice 
of parameters, the rescaled energy takes the form 


Wi,o(F; F) = tr [(F^F)^] - 2pi,o tr [F^F] + c^"^(F). (2.16) 


The next step of the proof is to show an affine relation between W^^^^ and Wi^. With Lemma 
|2.5[ we proceed by completing the square to get 

{R-,F) = tr [(F^F)2] - 2p^.^^ tr [F^F] + (F) 

= tr [(F^F)2] - 2p^.^^ tr [F^F] + + c^^}^^{F) 

= tr [(F^F - : (2-17) 

where cjj^,l^(F) := cj^,l^(F) — Inserting p = 1 and pc = 0, we obtain the special case 

Wi,o(F;F) =tr[(F^F-pi,ol)"] + c[^^(F) . (2.18) 


We can now reveal the connection between the minimization problem with parameters p > Pc > 0 
and the non-classical limit case (p,Pc) = (1,0). Note first that 


W^,^^{R;F)=tT 

\R^F-p^,^lf 

+ 41(7^)=tr 

(r^F-^ Pi,o l) 




V Pi,0 J 


+ 4%Jf) 


= tr 


1 ~ Pi,o 1 


+ 41 ( 7 ")- 


(2.19) 


In the last equation, we easily discover the trace term of equation (2.18) with one essential change: 
F was replaced by F^^^^ := A“^F. We now solve (2.18) for the trace term, and insert F^^^^ for 
the parameter F. This gives 


tr 


(F^F^,^^ - pi,ol)"l = lTi,o(i?;F^,^J - ■ 


( 2 . 20 ) 
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Next, we substitute the trace term in (2.19) by its expression in terms of Wi^o and finally obtain 

+41(F) 

= +cW^(F) , (2.21) 

with := This establishes the missing link ~ 

Wi^o(^; since > 0 and the final constant depends only on F. With this, the 

chain 

W^.^,(i?;F) ^ ~ (2.22) 

is now complete. All four energies give rise to the same energy-minimizing rotations in SO(n). I 


Once the optimal energy-minimizing rotations for Wi^o are available, the optimal rotations for 


3.4 


with general weights /r > /i^ 

This procedure is detailed in ! 

^ > /ic > 0, i.e., with strictly positive Cosserat couple modulus f^c, is completely governed by the 
case with zero Cosserat couple modulus fj,^ = 0 which is highly interesting in view of El! 


> 0 can be directly inferred by a substitution of F with 
In this sense, surprisingly, the non-classical case 


3 Optimal rotations for the Cosserat shear-stretch energy 


In this section, we compute explicit representations of optimal planar rotations for the Cosserat 
shear-stretch energy, i.e., we focus on dimension n = 2. The parameter reduction strategy 
in Lemma |2.2| allows us to concentrate our efforts towards the construction of explicit solutions 
to Problem |1.4| on two representative pairs of parameter values p. and pc- The classical regime 
is characterized by the limit case {p, Pc) = (1) 1) and the unique minimizer is given by the polar 
factor Rp(A) for any dimension n > 2, see Corollary 2.4 The non-classical case represented by 


{p, Pc) = (1)0) turns out to be much more interesting and we compute all global non-classical 
minimizers rpolar]^g(A) for n = 2. This is the main contribution of this section. Furthermore, 
we derive the associated reduced energy levels IT{®(^(F') and IT[®(5*(F) which are realized by the 
corresponding optimal Cosserat microrotations. Finally, we reconstruct the minimizing rotation 
angles for general values of p and pc from the classical and non-classical limit cases. 


3.1 Explicit solution for the classical parameter range: fic > ^ 0 


By Corollary 2.4 the polar factor Rp(F') is uniquely optimal for the classical parameter range in 


any dimension n > 2. Let us give an explicit representation for n = 2 in terms of Op € (—tt, tt]. In 


view of the parameter reduction, distilled in Lemma |2.2[ it suffices to compute the set of optimal 
rotation angles for the representative limit case {p,Pc) = (!> !)• 


Thus, to obtain an explicit representation of ap € (— tt, tt] which characterizes the polar factor 
Rp(A) in dimension n = 2, we consider 


argmin ITi^i(i?(a); F) = argmin 

a G [—tTjTt] a G [—7r,7r] 


COS a 

— sin a\ ^ 

{Fn 

Fi 2 

sin a 

cos a J 

[F 21 

F 22 



(3.1) 


Let us introduce the rotation J 0 / ^ SO(2). Its application to a vector u S 

corresponds to multiplication with the imaginary unit i € C. In what follows, the quantities 
tr [F] = Fii + F 22 and tr [JF] = —F 21 + F 12 play a particular role and we note the identity 

tr[F]^+ti[JFf = ||Ff-k2det[F] =tr[C/]^ . (3.2) 


The reduced energy IFi)i^(F) := min^ggo(„) Wip{R ; F) realized by the polar factor Rp(F) can be 
shown to be the euclidean distance of an arbitrary F in to SO(n). For n = 2, we obtain 
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Theorem 3.1 (Euclidean distance to planar rotations). Let F G GL“''(2), then 


SO(2)) = IIU - 111" = ||Ff - 2 ^||F|l" + 2det[F] + 2 . (3.3) 

The unique optimal rotation angle realizing this minimial energy level satisfies the equation 


/sin Op\ _ 1 /—tr [JF]\ 
l^cos ap) ~ tr [U] V tr [F] J 


In particular, we have otp{F) 


= arccos 



(3.4) 


Proof. See [H] [Appendix A2.1]. 


Corollary 3.2 (Explicit formula for Rp(F)). Let F G GL~''(2), then the polar factor Rp(E) has 
the explicit representation 


Rp{F) = R{ap) 


/cosOp — sinOpA _ 1 / tr[F] tr[JE]\ 

sin Op cos Op J ti[U] \-tr [ JF] tr [F] J 


(3.5) 


3.2 Symmetry of the first Cosserat deformation tensor 

The goal of this subsection is two-fold: first, we want to solve the equation skew(F'^F) = 0 for 
R G SO(2) which is equivalent to F^F G Sym(2). The unique polar factor Rp(F) is certainly a 
solution, but are there others? Second, we want to introduce an approach based on a rotation R 
relative to the polar factor Rp(F). This turns out to be essential to fully grasp the symmetry of 
the non-classical minimizers rpolar(F). This is the content of the next lemmaj^ 

Lemma 3.3 (Symmetry of the planar first Cosserat deformation tensor). Let F G GL'''(2) be 
given. The first Cosserat deformation tensor U{R) := R^F is symmetric if and only if 

F = ±Rp(F). (3.6) 


Proof. Let us first transform the equation into the orthogonal coordinate system induced by the 
principal directions of stretch. The orthogonal basis given by the eigendirections of U makes up 
the columns of a matrix Q G SO(2). As it turns out, it is natural to define a relative rotation in 
principal stretch coordinates which is given by 


F(/3) := Q^R{af Rp(F)Q = R{af Rp(F) . 


(3.7) 


Note that the rightmost equality holds only for SO(2), because it is commutative. We expand 
U := R^F = F^Rp(F)t7 = R^ Rp{F)QDQ'^ and exploit Q'^ sk.ew{X)Q = skew(Q^A'(5), i.e., the 
fact that skew is an isotropic tensor function. This gives 


skew([/) = 0 skew(F^ Rp(F)(5FQ^)(3 = 0 skew(Q^F^Rp(F)(5Fl) = 0 , (3.8) 

= :R 

where D := diag(cri, 02 ) ■= (^ i is the diagonalization of U. A simple computation in com- 

\ U 0-2/ 

ponents leads to the necessary and sufficient condition 


0 = ske'w{R{l3)D) 


( ° 

\^-i((Ti -kCT2)sm/3 


5(^1 


- (T 2 ) sin/3 

0 


sin/3 


tr[F] 




(3.9) 


We conclude that the necessary and sufficient condition for U G Sym(2) is sin(/3) = 0. Let us 
restrict /3 G (—tt, tt], then /3 = 0 V /3 = tt, i.e., R{P) = ±1. Substituting this into and solving 
for R{a) yields the claim. ■ 


®Cf. also |18l Eq. (2.12)] for the case n = 3. 









Remark 3.4 (Symmetry of strains vs. symmetry of stresses). Consider an energy W'^{U). Then, 

the symmetry of the Cauchy stress tensor (j{F) := 1 /(lei[F\ RDjjW^{U) F"'" is equivalent to the 

_ __ 

symmetry ofD-ijW^(U)U as was shown in m- We recall that the microstrain tensor U — I .2 

is symmetric if and only if R = ±Rp(i^). This symmetry does imply that the Cauchy stress 
tensor is symmetric. It is, however, possible that non-symmetric microstrains induce a 
symmetric Cauchy stress tensor. This may be unexpected, but it is the natural scenario for 
the case of non-classical optimal Cosserat rotations. A thorough discussion is given in m- 


3.3 The limit case {/j,, fic) = (1?0) for /j, f> Hc > 0 

We now approach the more interesting non-classical limit case {p, Pc) = (IjO) and compute the 


optimal rotations for W^^^^{R-,F). Note that, due to Lemma 
entire non-classical parameter range /r > /ic > 0. 


2.2 


this limit case represents the 


In the proof to Lemma 2.2 we have introduced Wi^q ~ Wi^o, i.e., a modified energy that gives rise 


to the same optimal rotations. In a similar spirit, inserting the particular values (/i, pc) = (1) 0) 
into (2.11) from the proof of Lemma 2.2 we can specialize to n = 2 as follows: 

WpoiRlF) = ltr[iR^F)^]-2tT[R^F]+^\\Ff + \\l2f 

^ ^tr[R'^F]‘^-2tT[R^F]+^\\Ff-det[F] + 2=:W{R;F) + c{F) . (3.10) 


= : W{R-,F) 


= ■■ &{F) 


This implies W ~ Wip, since both energies differ by a constant (with respect to R) 
c{F) = -det[i^]-b2 ^ ^tT[U]^ -2det[U]+2 . 


(3.11) 


The next step is to determine the critical rotations for F) by taking derivatives w.r.t. 

R S SO(2). Let us compute the necessary conditions. 

Theorem 3.5 (Characterization of the critical rotations for Wifi{R', F)). Let F G GL“''(2). A 
rotation R G SO(2) is a critical point for the energy Wifi{R; F) if and only if 

skew{R^F) = 0 , or tr =2 A tr [U] > 2 . 


Proof. Taking variations 6R = A ■ R, A G so{2), we arrive at the stationarity condition 

VAeso(2): (tr [R^F] - 2) ^) = 0 . (3.12) 


This equation holds good, if and only if either of the two factors on the left hand side vanishes. 

we have shown that skew(i?’^F) = 0, if and only if ii = ±Rp(F). Let us discuss 


3.3 


In Lemma 


the second possibility 


tr [R{afF] = 2 



(3.13) 


as an equation for a. The equation on the right hand side is easily obtained by a short computation 
in components. From the relation (■(;(«), w(F)) = cosa ||u(a)||||u>(F)|| = 2 it follows that the angle 
OL can only be solved for, if 


2 

Mm 


i/tr [F]Vtr [JFf 


(3.14) 


For tr [C/] > 2, due to the symmetry of the cosine, there exist two symmetric solutions for a 
(which may coincide). For 0 < tr [t/] <2 there is no solution. ■ 
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Let us enumerate the previously obtained necessary conditions for a critical point: 

1.) i?=-Rp(F), 2.) i? = +Rp(F), and 3.) tr [R^F] = 2 A tr [t/] > 2 . 


Note that for the two classical critical points =FRp(F), we have tr = =Ftr[[7]. Due to the 

particular expression of the energy W in terms of tr [ C7 ], we can insert the critical values of 
tr [R^F] into the defining equation ( |3.10[ ) for W to obtain the associated critical energy levels: 

= + 2tr[U], W^‘^\F) = \Y[Uf-2iv[U], and W^^\F) =-2 . 

Our first observation is that W^^\F) > W^‘^'>{F) for all F G GL“''(2). Further, if tr \U\ > 2, then 
the branch W^'^^F) exists and we have W^^\F) > W'^‘^\F) > IR*^^^(F), i.e., the non-classical 
branch realizes the global minimum when it exists. For 0 < tr [[/] < 2, the branch does 
not exist and the global minimum is realized by the classical branch W^‘^\ i.e., Rp(F) is uniquely 
optimal. Finally, for tr \U\ = 2, we have 

Theorem 3.6 (The formally reduced energy VF(®g'^(F)). Let F G GL“''(2), let the energies 

andW^^^ as above and W^'^^F) := (F)+ c(F), i = 1, 2, 3. Then, t/ie formally reduced 

energy 

Wr^(F):= min Wio(R;F):= min || sym(R^F - l)f (3.15) 

■ flGSO(2) ’ flGSO(2) 

is given by 


W^f{F) 


1F(2)(F) = tr [([/ - 1)2] = dist^(F, SO(2)) , 

TT(3)(f) = I ||Ff - det[F] ^ ^tr [R]^ - 2det[t/] , 


if tr [U] < 2 
if tr [R] > 2 . 


(3.16) 


Proof. It suffices to add the constant c(F) to the minimal energy levels for IF(R;F). 
Since W^‘^\F) corresponds to R = Rp(F) for which R^F = R is symmetric, we find that 
W^'^\F) = minflggo( 2 ) lRi,i(R; R) = dist^(F, SO(2)). Note that for tr [R] = 2, we have 

IF(2)(F) = VF(3)(i?). ■ 


It is well-known that any orthogonally invariant energy density IF(F) admits a representation in 
terms of the singular values of F, i.e., in the eigenvalues of R. Let us give this representation. 

Corollary 3.7 (Representation of VFi®o^(F) in the singular values of F). Let F G GL’*'(2) and 
denote its singular values by ai, i = 1,2. The representation of W[‘^q{F) in the singular values of 
F is given by 


W{':^{F) = Wlf{a^,a2) 


(cTi - 1)2 -I- (0-2 - if , if ai + <72 < 2. 
- a2f , if (7 i+(72> 2 . 


(3.17) 


Proof. We insert ||F||2 = ||R||2 = crj + cr| and det[F] = det[R] 
to see that both pieces of the energy coincide for cti + 0-2 = 2. 


(71(72 into (3.16). It is not hard 


Note that the previous formulae are independent of the enumeration of the singular values. 

3.3.1 Optimal relative rotations for /j, — 1 and /Xc = 0 

Our next goal is to compute explicit representations of the rotations rpolar^g(F) which realize the 
minimal energy level IF^^^(F) in the non-classical limit case (/i,/T c) = (IjO). This is the content 
of the next theorem for which we now prepare the stage with the following 

Lemma 3.8. Let D = diag((7i, ( 72 ) > 0, i.e, a diagonal matrix with strictly positive diagonal 
entries. Then, assuming tr [D] > 2, the equation tr [R(/3) R] = 2 has the following solutions 

13^ = ± arccos ^ [-tt.tt] . (3.18) 

For tr [D] < 2, there exists no solution, but we can define /3 = f3^ := 0 by continuous extension. 
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/3to(tr[C/]) 



Figure 3.1: Plot of the two optimal relative rotation angles (3^^ = iarccos for the non- 

classical limit case (^,/rc) = (1)0)- Note the pitchfork bifurcation in tr [17] = pi_o = 2. For 
0 < tr [U] <2, the polar angle Ofp is uniquely optimal and the relative rotation angle j3 vanishes 
identically. 


Proof. We compute 


/cos P — sin /3\ /(Ti 0 \ 
\^sin/3 cos/3 / \ 0 cr 2 / 


(cti -I- 0 - 2 ) cos /3 = tr [D] cos /3 . 


(3.19) 


Since tr [U] > 0, we may divide to obtain the relation cos/3 = 2/tr[i7]. This is solvable if 
and only if 2/tr[i7] < 1 which is equivalent to tr [13] > 2. There are two symmetric solutions 
= ± arccos(2/tr [73]). Since both vanish for tr [73] = 2, we can continously extend /3 = /3^ := 0 
for tr [73] <2. ■ 


Our Figure [3T] shows a plot of the optimal relative rotation angle ^(tr [[/]). In the classical parame¬ 
ter range 0 < tr [17] < 2, Q;p(7^) is uniquely optimal and P vanishes identically. In tr [17] = 2, a clas¬ 
sical pitchfork bifurcation occurs. In particular, due to tr [17(12)] = tr [I 2 ] = 2, the identity matrix 
is a bifurcation point of P^{F). Further, we note that the branches /3^(tr [17]) = ± arccos(2/tr [17]) 
are not differentiable at tr [17] = 2. This has implications on the interaction of the Cosserat 
shear-stretch energy with the Cosserat curvature energy IFcurv 

Theorem 3.9 (Optimal non-classical microrotation angles afo)- ^ ^ GL'''(2) and consider 
(/i, /Xc) = (1, 0). The optimal rotation angles for Wi^o given by 


Q!p(F) = arccos(^) , if tr [17] < 2 

ap{F) ± arccos (tTpj) , if tr [17] > 2 . 


(3.20) 


3.6 


3.5 


which shows that Rp(7^) realizes Fur- 

(see also the proof) the branch uniquely corresponds to Rp(F). In 


Proof. The first statement follows from Theorem 
ther, by Theorem 
other words 


tr[[/]<2 ^ a^,^^(7^) = ap(F) rpolar(F) = Rp(F) . (3.21) 

Let us now assume tr [U] > 2. In this case, by Theorem |3.5[ globally energy minimizing rotations R 
realize Thus, a G (—TT, tt] is a solution of tr [ R{a)^ ^F] = 2 for given F G GL+(2). Consider 

again the relative rotation (cf. the proof of Lemma |3.3[ ) given by 

R{P) ■.= R{af Rp{F) . (3.22) 

The uniqueness of Rp(F) implies a one-to-one correspondence between R{a) and R{P). In terms 
of rotation angles, we find that 

a = ap — P , (3.23) 
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where ap denotes the rotation angle of the polar rotation Rp(_F), i.e., R{ap) = Rp(J^). After a 
transformation into the coordinate system given by the principal directions of stretch (i.e., given 
by the eigendirections of C/), we obtain 


2 = tr [R^F] = tr R{I3)D . (3.24) 

Applying Lemma|3.8| we find that there are two energy-minimizing relative rotation angles 


?^ = ± arccos 


tr [D] 


= ± arccos 


tr[C/] 


, for tr [U]> 2 . 


We can now solve (3.231 for the corresponding microrotation angles a^Q which gives 

- 2 

'^10 = Oip ~ P = ttp T arccos 


tr [C7] 


(3.25) 


(3.26) 


The second equality is just another application of Lemma B 


3.4 General values for ^ and Hc 


The reduction for /i and /ic in Lemma |2.2| asserts that the optimal rotations for arbitrary values 
of ^ > 0 and fic>0 can be reconstructed from the limit cases (^, fic) = (1,1) and {fi, /ic) = (1,0). 
We now detail this procedure which essentially exploits Definition |2.1[ 

Note first that the rescaled deformation gradient := induces a rescaled stretch tensor 




(3.27) 


The right polar decomposition takes the form = Rp(T’p_pJ {7p,p^. From Rp(Tp^pJ = 

Fp_p^t/“p^ follows the scaling invariance Rp(F)i,/ic) = Rp(A). For the non-classical parameter 
range /i > /ic > 0, the quantity 


tr 


U, 




Pp.Pc 


plays an essential role. This leads us to 


tr 


U, 




^ 2 — pi Q 


tr 


Pi,o 

Pp.Po 


U 


> Pi.o 


tr [U] > pp,p 


(3.28) 


(3.29) 


In particular, this implies that the bifurcation in tr \U] allowing for non-classical optimal planar 
rotations is characterized by the singular radius Pp,p,, := 

Theorem 3.10. Let F G GL~'’(2). For /ic > /r > 0 the optimal microrotation angle is given by 


— rrp(A) — arccos 


For fJ. > Pc ^ 0, the two optimal rotation angles are given by 


tr[f] 

tr[t/] 




p(F) = arccos , if tr [U] < Pp, 

ap{F) =F arccos , if tr [17] > Pp 


Me 
— PM^Mc 


(3.30) 


(3.31) 


Proof. By Corollary 2.4 Rp(F) is uniquely optimal for the classical parameter range /ic > /i > 0. 
The associated rotation angle ap{F) is immediately obtained from the explicit formula for the polar 


factor given in Corollary |3.2[ Let^us now discuss the more delicate non-classical parameter regime 
M > Me > 0. Here, the rescaling J^p,p^ plays a decisive role. First, the condition 0 < tr [17] < Pp,p,, 

is equivalent to 0 < tr 17p_p^ < 2. In this case, the polar factor Rp(Ap,p^) = Rp(A) is uniquely 
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optimal. For the parameter domain tr [U] > which is equivalent to tr > 2, however, 

we obtain the optimal relative rotation angles 


PtAP) = = ± arccos 


Pi.o 


tr 


Uu, 


= ± arccos 


Pm.M e 

tr[U] 


(3.32) 


This gives ai^o(^M.Mc) = o^piAA “ AoiAA = “ AoiAA- 


4 Optimal rotations for planar simple shear 


We now apply our previous optimality results to simple shear deformations. Previously, in |23| . 
Neff and Miinch contributed the optimal planar rotations for simple shear. A simple shear of 
amount 7 S M is a homogeneous linear deformation represented by a matrix of the form 

F, := (; ]) . (4.1) 

In this section we derive the energy-minimizing rotation angles := for simple 

shear. 

Let us shortly consider the classical limit case {p.,Pc) = (1.1) which represents the entire classical 
parameter range p-c > p > 0. Essentially due to Theorem o the polar rotation Rp(F.y) is 
then uniquely optimal, see Corollary |2.4[ Thus, we proceed with the non-classical limit case 
(p. Pc) = ( 1 . 0 ) which represents the entire non-classical parameter range p > pc > 0 , as we have 
seen in Lemma 12.21 

Let us collect some properties of simple shear F^. We have ||Fly||^ = 24 - 7 ^ and det[E.y] = 1, i.e., 
simple shear is volume preserving for any amount 7 . This allows us to compute 

tr [U^] = + 2det[E^] = > 2 . (4.2) 

Thus, the reduced energy always satisfies IFi®o*(E. 7 ) = for simple shear F^, i.e., the 

non-classical branch is always optimal. 

Corollary 4.1 (Optimal non-classical Cosserat rotations for simple shear). Let (p,Pc) = (1,0) 
and let F^ G GL^(2) be a simple shear of amount 7 S K. Then, 

7^0 rpolar)^o(F..y) ^ Rp(FV) . (4.3) 


Proof. First, tr [?7.y] > 2 for all 7 S K and with Theorem 3.5 the optimal relative rotation angle 
(3 G [—TT, tt] satisfies 

For 7 ^ 0, it is easy to see that tr [UA > 4. Since arccos(2/a;) is strictly increasing for a; > 2, we 
finally conclude: 

2 2 

0 < \/3{F.y)\ = arccos( ^ rrri ) “ arccos(- 


^ tr [U.y 


'\/4T 


T 


()• 


Remark 4.2 (Symmetry of the first Cosserat deformation tensor U in simple shear). A simple 
shear Fj by a non-zero amount 7 0 automatically generates an optimal microrotational response 

rpolar^(i4y) which deviates from the continuum rotation Rp(F'). This implies that the associated 
first Cosserat deformation tensor ^{F.^) rpolar)*^Q(Fly)^E.y is not symmetric for any 7 ^ 0; 
cf. Lemma\3^ 


Remark 

„± 


4.3 (Consistency with |23jL It is not hard to show that the explicit minimizers 
rpolarf Q(Fly) for the optimal Cosserat rotations previously obtained do exactly match those com¬ 
puted in 12, Equation (3.24-)]. We have found the following identity to he helpful for the 

verification: arctan(^) = sign( 7 ) arccos( ^ ) = sign( 7 ) arccos( ^.^i-^ j ) . 
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Figure 4.1: Plot of the critical energy levels i = 1,2,3, of Wi_o for a simple shear 

of amount 7 . Note that The critical energy levels are realized by 

1 F(i) = — = +Rp(F.y) and IF^^^ = rpolari^g(F^), respectively. 


-10 -5 0 5 10 
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Figure 4.2: Plot of the optimal microrotation angles a!i’^g( 7 ) for Wi^o and simple shear F^ of 
amount 7 S M. The negative optimal branch a!j”g( 7 ) [dashed green curve] exactly eliminates 
the angle Q:p(F) and vanishes identically. The positive branch Q!]j^g( 7 ) [continuous blue curve] 
corresponds to a rotation by 2 ap( 7 ) [dot-dashed red curve]. Note the symmetry w.r.t. to the 
continuum rotation angle Q:p( 7 ) 


4.1 Simple glide and cancellation of the polar factor 

For the case of simple shear, one of the optimal Cosserat rotations for the shear-stretch energy 
Wi.o exactly cancels the polar factor. More precisely, one of the two rotations rpolar^(Fy) is 
the identity element I 2 S SO(2), while the other solution is given by Rp(F!y)^, i.e., a rotation by 


2ap(Fg,), see also Figure 4.2 It would be quite intriguing if this “gliding” behavior were specific 


to simple shear, but, as it turns out, it is possible to construct other examples showing the same 
behavior. 

To see this, note first that 


(Cor. 

ttn — — arccos 


tr[F] 

tr[C/] 


and 1/3] arccos 

' Vtr[C/] 


(4.5) 
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The condition that one of the relative rotations cancels the continuum rotation is given by 


l/5| = |ap| 


tr[F] 

tr [17] 


2 

Wi’ 


(4.6) 


whenever tr[U] > 2. Hence the set of matrices for which the polar rotation is canceled by a 
minimizing relative rotation is given by {F G GL'''(2) | tr [F] = 2 A tr [C7] > 2}. This set is 
nonempty, because it contains the family of simple shears Fj. In order to see that this set also 
contains homogeneous deformations which are not simple shears, we consider the matrix 



and set 

Fj,k := Fy + k1 * and '■= ^Fy^^ . 

Note that det)^^,^] = 1 — implies that G GL^(2) for 0 < |k| < 1. Further, 
Vy, K e M : tr [F^.k] = tr [F.y] and tr [JF.y_„] = tr [JF.y] , 

which implies 

tr [Uy^f^f = tr [F.y,„]^ + tr [J(F.y_«,)]^ = tr [F.y]^ + tr [JFyf = tr [17^,]^ . 


(4.7) 

(4.8) 


(4.9) 


(4.10) 


Hence, both quantities tr [LTy^fj] and tr [F^^k] are independent of n and condition (4.6) is automat¬ 
ically satisfied for all admissible F^y^^, i-e., for 0 < |/t| < 1 . 


We conclude that, given any simple shear Fy of amount 7 S K, there is a one parameter family 
Fly,K, 0 < |k| < 1 of matrices that are not simple shears for which one of the optimal relative rota¬ 
tions R exactly cancels the continuum rotation Rp(F). The interesting “glide behavior” observed 
in Section]^ is not specific to simple shear. 


5 Conclusion 


In Section we have seen that it is sufficient to construct energy-minimizing rotations for the 
classical limit-case (/r,/ic) = (IjO) and the non-classical limit case (/r,^c) = ( 1 ) 1 )? respectively. 
For the classical parameter range p.c > M > 0) the unique minimizing rotation for F) 

is given by Rp(F), in any dimension n. For ^ = /Tc, the reduced Cosserat shear energy can be 
formally reduced to 

WltiF) = ; F) = ■ (5-1) 

Hence, setting the Cosserat curvature coefficient = 0, one can expect the full quadratic Cosserat 
model to behave essentially like a classical Biot model, see, e.g., the introduction to |I 8 j . 


However, a fundamental motivation to use extended continuum models such as Cosserat models, 
is to generate solutions showing non-classical effects. For the quadratic Cosserat model (without 
curvature), this is the case if there is a deviation R ^ Rp(F), since the model (formally) reduces 
to the well-known Biot energy otherwise. In Section we have shown that this is only to be 
expected for the non-classical parameter range /r > /ic > 0. If non-classical solutions should be 
generated already in the identity I 2 , then we even have to require /ic = 0 , since > 2 otherwise. 
The existence of the presented non-classical minimizers rpolar(F) is highly interesti ng. In strong 


contrast, if we replace the non-symmetric strain tensor t/ — 1 by log U in Problem 1.2 which is 


natural in view of the Cartan decomposition of 0 ((?t,), one can show that no non-classical solutions 
exist for arbitrary /i > 0 and /ic > 0 


argmin |/i ||symlog(F^F)||^ + pc ||skewlog(F^F)||^| = {Rp(F)} . (5.2) 

Resoin) ^ ^ 


For a proof and a deep discussion of the nature and properties of logarithmic strain measures, 

see [IIlElEl]. 

In our introduction, we have stated that the solution to Problem |1.2| for n > 3 is currently out of 
reach. However, we have successfully computed non-classical critical Cosserat microrotations for 
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n = 3 using a parametrisation by unit quaternions and computational algebra. Further, we have 
managed to select the energy minimal branches, experimentally. An extensive numerical validation 
shows, moreover, that our candidates are very likely the global minimizers. The mechanisms 
discovered for the case n = 2 in the present work do carry over to the case n = 3 quite literally up 
to the determination of the microrotation axis. This is the content of a forthcoming second part 
of this paper jS]. For dimensions n > 3, the weighted Problem 1.2 is, to the best of our knowledge, 
still completely open. It seems to us, however, reasonable to guess that a transformation into the 
principal directions of stretch, i.e., the eigendirections of U, is a good plan of attack. 


Remark 5.1 (Final Conclusion). To ascertain the complete absence of a non-classical response 
within a geometrically nonlinear quadratic Cosserat-micropolar shear-stretch energy, one must 
choose a classical parameter set, i.e., /ic > /r > 0. 
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A Appendix 

A.l Some planar matrix identities for U and F 

Let n = 2. Applying the trace to both sides of the Cayley-Hamilton equation and exploiting 
linearity, we obtain 

tr - tr [X] X + det[A] I 2 ] = tr [0] ^ tr [X^] = tr [X]^ - 2 det[X] . (A.l) 

This leads us to the following identity. 

Lemma A.l. Let F G GL''’(2) and U := V F'^F G PSym(2). Then, 

tr [U] = ^J\\Uf + 2det[U] = ^J\\Ff + 2det[F] . (A.2) 

Proof. Note first that tr [t/^] = (t/, U) = \\U\\'^ and that tr[t/] > 0. The expression in terms 
ofF e GL+(2) is implied by det[F] = det[Rp(F)^F] = det[C/] and IJF]]^ = ||Rp(F)'^Ff = \\uf.m 
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